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ABSTRACT 

The recent Fermi detection of the globular cluster (GC) 47 Tucanae high- 
lighted the importance of modeling collective gamma-ray emission of millisecond 
pulsars (MSPs) in GCs. Steady flux from such populations is also expected in the 
very high energy (VHE) domain covered by ground-based Cherenkov telescopes. 
We present pulsed curvature radiation (CR) as well as unpulsed inverse Compton 
(IC) calculations for an ensemble of MSPs in the GCs 47 Tucanae and Terzan 5. 
We demonstrate that the CR from these GCs should be easily detectable for 
Fermi , while constraints on the total number of MSPs and the nebular B-field 
may be derived using the IC flux components. 

Subject headings: gamma rays: theory — globular clusters: individual (47 Tu- 
canae, Terzan 5) — pulsars: general — radiation mechanisms: non-thermal 


1. Introduction 

The launch of Fermi Gamma-ray Space Telescope (Fermi) on 11 June 2008 heralded 
a new era for high-energy (HE) Gamma-ray Astronomy. The Large Area Telescope (LAT) 
aboard Fermi (Atwood et al. 2009) is ~ 25 times more sensitive than its predecessor, the 
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Energetic Gamma Ray Experiment Telescope (EGRET). Recently, Fermi anounced the de- 
tection of the globular cluster (GC) 47 Tucanae as a bright point source (Abdo et al. 2009). 
This emission may plausibly be due to the large number of millisecond pulsars (MSPs) 
residing in this cluster (e.g. Harding, Usov, & Muslimov 2005, hereafter HUM05). 

GCs are furthermore possible sources of very high energy (VHE) gamma rays (Bednarek 
& Sitarek 2007, hereafter BS07), and are / will be important targets for ground-based 
Cherenkov telescopes such as the High Energy Stereoscopic System ( H.E.S.S . - Aharonian 
et al. 2006) as well as future telescopes such as the Cherenkov Telescope Array ( CTA ) l and 
the Advanced Gamma-ray Imaging System (AGIS) 2 . 

Venter & de Jager (2008a) modeled the collective curvature radiation (CR) expected 
from 47 Tucanae, demonstrating that this pulsed component should be detected by Fermi , 
depending on the number of visible gamma-ray pulsars in the GC MSP population, N vis (see 
§2). Venter & de Jager (2008b) extended these results by modeling the unpulsed inverse 
Compton (IC) and synchrotron radiation (SR) flux components expected due to the interac- 
tion of HE electrons, ejected from MSP magnetospheres, with the ambient cosmic microwave 
background (CMB) and bright starlight photons, in the presence of a nebular magnetic field 
B. 


In this Letter, we give updated CR spectra (§3, see Venter & de Jager (2009) for details), 
in addition to further IC calculations (§4), for 47 Tucanae and also Terzan 5, another GC 
containing many MSPs (§2). The roles of Fermi, H.E.S.S., and CTA for deriving model 
constraints are discussed in §5. 


2. Globular Cluster Millisecond Pulsars 

GCs are very old galactic sub- structures, with evolved stellar composition. As a con- 
sequence, they are likely to harbor a larger than usual density of end products of stellar 
evolution - compact objects. This is supported by observations of X-ray binary systems 
and MSP populations in GCs. Since the discovery of the first GC MSP in M28 (Lyne et al. 
1987), a total of 140 pulsars have been discovered in 26 GCs (Freire 2009), with 90% having 
periods P < 20 ms (Ransom 2008). 

Among the GCs, 47 Tucanae and Terzan 5 stand out by the number of radio-detected 

1 www. cta-observatory.org 

2 http: / /gammal. astro. ucla.edu/agis/index.php/Main_Page 
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MSPs they host: 23 and 33 respectively (Ransom 2008). However, these GCs are other- 
wise remarkable for different reasons. 47 Tucanae is one of the most massive GCs, after 
lo Centauri, while the core density of Terzan 5 gives it the highest expected rate of binary 
interactions (Pooley & Hut 2006). Basic properties for these two objects are summarized in 
Table 1. 

Besides pulsed emission, detected already in radio, X-rays, and gamma rays, MSPs are 
also expected to produce steady low- flux VHE emission. Given the sensitivity of current 
instruments, it is likely that only a population of such objects would produce detectable 
VHE fluxes (§4). The extension of the cores of the GCs is furthermore clearly below the 
angular resolution of current HE and VHE instruments. In this Letter we are therefore 
not concerned about individual MSP fluxes, but attempt to constrain a particular class of 
polar cap (PC) pulsar models using the total flux from a GC MSP ensemble. This approach 
reduces uncertainties in the CR and IC fluxes due to unknown pulsar geometries, as the 
relative uncertainties in ensemble-averaged fluxes, which scale as the reciprocal of the square 
root of the number of MSPs (Venter & de Jager 2008a), are typically one order of magnitude 
smaller (for ~ 100 MSPs) than for the case of single MSP ‘geometry- averaged’ fluxes. 

An important issue in the context of this work is the uncertainty on the size of the 
MSP population in each system, beyond the radio-detected ones. We dinstinguish between 
two numbers. We define the number of visible gamma-ray pulsars (iV vis ) as those MSPs 
whose CR is (mostly) beamed toward earth, so that observers will in principle be able to 
measure pulsed emission from these (although it might not always be possible to distinguish 
individual members due to experimental limitations). In contrast, Wot represents the total 
number of MSPs in the GC, whether visible in pulsed gamma rays or not. 

Mostly off-axis pulsars might still contribute to the total CR flux, as observers may be 
clipping the faint edges of gamma-ray beams from these pulsars. We therefore include them 
in N v i S , along with the on-beam MSPs. HUM05 noted that gamma-ray beams of MSPs are 
probably wider than their radio beams, so N vis is probably larger when applied to gamma-ray 
pulsars (as done here) than when applied to radio ones, possibly up to N vis ~ Wot- The IC 
flux however depends on Wot, as all pulsars, including ones with off-beam geometries, are 
expected to contribute VHE leptons (mostly electrons, as most MSPs are believed to be pair 
production starved; see HUM05) which upscatter soft photons to gamma-ray energies. 

The isolated location of 47 Tucanae in the Southern sky away from the Galactic Plane 
made it a prime target for GC studies. In radio, McConnell et al. (2004) estimated from 
the unresolved emission a maximum of 30 detectable MSPs. In X-ray, Chandra observations 
revealed a large number of unidentified sources in the core of 47 Tucanae (Heinke et al. 
2005), from which followed an upper limit of about 60 MSPs, based on the X-ray properties 
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of radio-confirmed MSPs in the GC. 

For Terzan 5, the available observations are not so constraining. Again in the radio, 
Fruchter & Goss (2000) estimated a range of 60 — 200 MSPs in the core of the GC, from Very 
Large Array (VLA) and Australia Telescope Compact Array (ATCA) data. No equivalent 
result has yet been derived in X-rays. The analysis of Chandra data by Heinke et al. (2006) 
yielded X-ray sources, but with too high a sensitivity limit to constrain the MSP population. 

Finally, for both systems, GC evolution models (e.g. Ivanova, Heinke, & Rasio 2008) 
suggest MSP numbers below 200, of which a fraction would be detectable. 

3. Pulsed Gamma-ray Flux 

As previously (Venter & de Jager 2008a), we use an isolated pulsar PC model (Venter 
& de Jager 2005), including the effect of General Relativistic (GR) frame dragging (e.g. 
Muslimov & Harding 1997; Harding & Muslimov 1998). In accordance with the expectation 
that most MSPs are inefficiently screened by CR and IC scattering pairs (HUM05), we use 
an unscreened acceleration potential. We furthermore use the same population of MSPs: P- 
values (period-derivatives) from Bogdanov et al. (2006), which imply an average spin-down 
luminosity of E mt = 2.2 x 10 34 ergs -1 , as well as canonical values for equation of state (EOS) 
parameters (moment of inertia Ins = 0.4A7 R 2 . stellar radius R = 10 6 cm, and pulsar mass 
M = IAMq). As previously, we add the CR spectra of iV vis = 100 MSPs, randomly selected 
from our population of 13 members (each with a random pulsar geometry, i.e. magnetic 
inclination and observer angle, including off-beam cases) iteratively a million times, and so 
obtain an average pulsed spectrum, with 2<j-bands indicating the uncertainty due to varying 
P and P- values, as well as the unknown pulsar geometries. 

Venter & de Jager (2008a) used a delta function approximation for the CR spectrum 
radiated per primary, but the resulting spectral shape did not correspond well to results of 
similar studies (e.g. Frackowiak & Rudak 2005), as pointed out by Dyks (2008, personal com- 
munication). We therefore developed a corrective technique involving the full CR spectrum 
per primary (Venter & de Jager 2009), and applied it here. 

Assuming that the MSPs in Terzan 5 have similar basic pulsar parameters than those in 
47 Tucanae, we scaled the cumulative CR spectrum of the latter to distances of d = 5.5 kpc 
and 8.7 kpc. The resulting updated CR spectra for 47 Tucanae and Terzan 5 are shown in 
panels (a) and (b) of Figure 1. 
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4. Unpulsed Gamma-ray Flux 

We calculate the ejection spectrum of electrons leaving each MSP by binning the num- 
ber of primary electrons ejected per unit time according to their residual energies at the 
light cylinder. We add the spectra of N tot = 100 MSPs (with random inclination angles) 
and iterate the Monte Carlo procedure a million times. Prom this we obtain an average 
cumulative ejection particle spectrum and 2cr-bands (Venter 2008). 

We divide the region where unpulsed radiation is generated into two zones: Z 0, reaching 
from r = 0 to r = r c , with r c the core radius, and Z 1, reaching from r = r c to r = rh m , 
with r hm the half mass radius (see Table 1). For 47 Tucanae, we used energy densities of 
u rad ~ 3 000 eV/cm 3 for Z0 , and ~ 100 eV/cm 3 for Z 1 for the starlight component, 
assuming a temperature of T = 4 500 K (more details may be found in Venter & de Jager 
2008b). These become ~ 1 000 eV/cm 5 * 3 and u 3 ad ~ 40 eV/cm 3 for Terzan 5 (note that 
the energy densities, and therefore the IC ‘bumps’ associated with the upscattering of stellar 
photons, are linearly dependent on the assumed total stellar luminosity - see Table 1). For 
the CMB component, we use Ucmb ~ 0.27 eV/cm 3 (for T = 2.76 K). 

We next calculate the steady-state injection electron spectrum for each zone by multi- 
plying the cumulative ejection spectrum by an effective time-scale which incorporates both 
radiation losses and particle escape. We modeled the latter (i.e. trapping of particles by the 
nebular field) assuming Bohm diffusion. The injection spectra were calculated for nebular 
fields of B = 1/xG and B — 10/:tG (BS07). Using these results, we obtained the expected 
unpulsed IC flux for each GC (see Figure 1). While the lower-energy SR spectra are not dis- 
cussed here, we do include SR losses in the calculation of the steady-state injection particle 
spectrum (see Venter & de Jager 2008b). 

5. Discussion and Conclusions 

We have modeled the HE and VHE components of CR and IC flux expected from a 

population of MSPs in the GCs 47 Tucanae and Terzan 5, for 100 members (Figure 1). For 
illustration, we assume that all MSPs are visible (i.e. N v is = N tot = 100). However, the CR 
spectra are linearly dependent on N vis , while the IC spectra are linearly dependent on Wot- 

Fermi has recently released a list of bright sources (> 10o\ Abdo et al. 2009), including 
preliminary fluxes of 47 Tucanae of (5.1T2.0) x 10 -8 cm -2 s -1 and (5.6±1.0) x 10~ 9 cm'V 1 
for energy bands spanning 100 MeV - 1 GeV and 1 GeV - 100 GeV respectively. Our average 

CR spectrum is overpredicting these fluxes by factors of ~ 2 and ~ 7. However, the pulsed 
spectrum may quite reasonably be scaled to N vis = 50 (as suggested by observations). Also, 
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extrapolating from observed trends in single MSP spectra, one may speculate that it is 
possible that the cut-off energy may be smaller due to uncertainties in EOS parameters, 
the electric potential, magnetospheric structure, and pulsar geometry. For instance, taking 
N vis = 50, and halving all CR spectral energies (i.e. shifting the CR spectrum to the left by 
a factor of 2), the model predictions improve to ~ 0.5 and ~ 1.7 times the Fermi flux in 
these bands. 

Our ensemble-averaged CR spectrum, resulting from 10 8 single MSP spectra (i.e. 10 6 
100-member cumulative spectra) with typical spectral cut-offs around a few GeV (generally 
depending on the individual MSP’s P , P, EOS, and geometry), exhibits a broad peak around 
~ 3 GeV. This average spectral cut-off energy follows directly from the assumed acceleration 
potential and B-field geometry of the PC model under consideration, as well as the fact 
that we simulate a random selection from the P — P parameter space for GC MSPs (similar 
to HUM05), and with different geometries. Given the fact that our CR spectrum seems 
to be reasonably successful at reproducing the preliminary Fermi data, we expect that the 
predicted cut-off energy must be close to the real one. This implies that the average model 
acceleration potential must be of the right order of magnitude, assuming dipolar B-fields. A 
more detailed study involving further model and parameter space constraints will have to 
await further Fermi spectral results. 

We indicate EGRET upper limits (ULs) from Fierro et al. (1995) on Figure la, and from 
Michelson et al. (1994) on Figure lb, converted assuming an E~ 2 spectrum. Interestingly, 
the 100 MeV UL is close to the Fermi sensitivity for 47 Tucanae. The EGRET UL for 
Terzan 5 is not very constraining, and implies N vis < 363 for d = 5.5 kpc, and V vis < 907 
for d = 8.7 kpc (using the average integral CR spectrum). 

As a reference, we have also included a CR prediction from HUM05, scaled to N vis = 100, 
in Figure 1. HUM05 noted that the collective radiation from MSPs in GCs may be visible for 
Fermi , and proceeded to apply their unscreened GR frame-dragging PC model to the MSPs 
in 47 Tucanae. They estimated the CR flux above 100 MeV from 47 Tucanae by integrating 
along a single magnetospheric B-field line, normalising to a (conservative) solid angle of 1 sr 
per PC, and scaling by the standard Goldreich- Julian PC current. This approach, which was 
followed to circumvent full 3D modeling, may be regarded as yielding a theoretical upper 
limit, given that different pulsar geometries (especially off-beam ones) have not expicitly 
been taken into account as done here, thus significantly exceeding our results at 100 MeV. 

Predictions from BS07 focused on the possibility that HE leptons, which may be ac- 
celerated inside MSP magnetospheres and / or reaccelerated at colliding pulsar and stellar 
wind shocks, gradually diffuse through the GC and upscatter CMB and soft stellar photons 
to produce an unpulsed component in the GeV-TeV energy band. BS07 assumed that 1% 
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of the spin-down power (with the population average taken as E rot ~ 1.2 x 10 35 ergs -1 ) is 
injected in relativistic leptons. In the absence of detailed model predictions, they assumed 
power-law injection spectra (typical for shock acceleration) with various low and high-energy 
cut-offs and spectral indices, and normalized these using the assumed average MSP spin- 
down power and conversion efficiency. Representative resulting IC spectra, linearly scaled 
to our population-average spin-down luminosity of E TOt = 2.2 x 10 34 ergs -1 , are shown in 
Figure 1. Our IC predictions, obtained by conservatively assuming no reacceleration of elec- 
trons, roughly agree with those of BS07 around a few TeV, but deviate at other energies due 
to different injection spectra, as well as slightly different assumptions regarding the energy 
density of stellar photons in the GCs (compare our two-zone approach with BS07’s energy 
density profile given in their Figure 1). 

Constraints on N tot as well as B may be derived from our model, using the predicted 
IC flux and assuming non-detection above 1 TeV at the level of 1% and 0.1% of the Crab 
flux, as achievable by H.E.S.S. and CTA, respectively. As shown in Figure 2, regions of the 
B — iVtot parameter space could be excluded from VHE ULs, e.g. choosing 5=10 for 
47 Tucanae, the H.E.S.S. UL constrains the total number of MSPs to < 33 (and < 3 for 
CTA). Terzan 5 is less constrained (these limits become < 141 and < 14 for d = 5.5 kpc, 
and < 353 and < 35 at d = 8.7 kpc). For small 5, particle escape is greater (and IC flux is 
lower), so that the constraints are less severe. This is also true for large B , where higher SR 
losses lead to inhibited IC flux. Furthermore, taking reacceleration of VHE electrons into 
account would increase SR and IC fluxes, implying more stringent constraints on N tot and 
B in the case of a non-detection. 

The B — N tot constraints are however dependent on the choice of diffusion coefficient 
(a larger coefficient would result in enhanced particle escape and thus lower IC flux). It 
is thus possible to invert the argument and obtain constraints on the diffusion coefficient 
in the absence of reacceleration, using future VHE limits and assuming reasonable values 
for Mot and B. To this end, a combination of Fermi data and Monte Carlo modeling 
may help constrain N V - 1S and N vis /N tot , and thus N tot . Also, measurements of the predicted 
diffuse ultraviolet SR (Venter & de Jager 2008b) and unpulsed diffuse gamma-ray IC spectra 
may further constrain the lepton injection spectrum as well as B. Lastly, if B is large 
enough (> 20 /j,G) so that diffusion losses are unimportant relative to SR losses, particle 
reacceleration in the GC may be studied without having the concern of the latter being 
masked by energy-dependent diffusion. 
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sions. This research was supported by an appointment to the NASA Postdoctoral Program 
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Table 1. Selected Globular Cluster Characteristics 



47 Tucanae 

Reference 

Terzan 5 

Reference 

Distance d (kpc) 

4.0 

1 

5.5 - 8.7 

2,3 

Mass (Mq) 

1.0 x 10 6 

4 

3.5 x 10 5 

5 

Core radius r c (arcmin) 

0.40 

6 

0.18 

6 

Half mass radius (arcmin) 

2.79 

6 

0.83 

6 

Detected radio MSPs 

23 

7 

33 

7 

Total stellar luminosity (L©) 

7.5 x 10 5 

8 

1.5 x 10 5 

8 


References. — (1) McLaughlin et al. (2006); (2) Ortolani et al. (2007); (3) Cohn 
et al. (2002); (4) McLaughlin & Van der Marel (2005); (5) Ivanova, Heinke, & Rasio 
(2008); (6) Castellani (2009), extended from Harris (1996); (7) Freire (2009); (8) 
BS07. 
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log 10 [Energy /TeV] 

Fig. 1. — Differential CR and IC spectra for 47 Tucanae (la) and Terzan 5 (lb). The top 
IC spectra are for a nebular field B = 10 /iG, while the bottom ones are for B — 1 pG. Dark 
gray lines indicate Fermi (Atwood et al. 2009) and H.E.S.S. (Hinton 2004) sensitivities 
(appropriate for each GC). Arrows represent EGRET upper limits. In (la), we indicate 
a scaled prediction from HUM05, as well as 2a-bands (light gray bands) due to unknown 
pulsar geometry (the bottom CR one stops where 2 a exceeds the mean value). For both 
panels, we show predictions from BS07 (scaled to an average spin-down luminosity of E mt — 
2.2 x 10 34 ergs -1 ). In (lb), solid lines are for d = 5.5 kpc and dot-dashed ones for d — 8.7 kpc 
(2a-bands not shown). 



- 12 - 



5 10 15 20 25 30 35 40 45 50 

Nebular B-field (jaG) 

Fig. 2. — Constraints on the total number of MSPs N tot and nebular magnetic field B. 
Panel (2a) is for 47 Tucanae, (2b) for Terzan 5. The dark gray indicate regions excluded 
by a H.E.S.S. non-detection, while light (plus dark) gray regions are excluded by a CTA 
non-detection (see text for details). In panel (2b), the solid lines are for d = 5.5 kpc, and 
dashed lines for d = 8.7 kpc. The thick horizontal lines indicate the number of currently 
detected radio MSPs (see Table 1). 


